Phase fluctuations of s-wave superconductors on a lattice 
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Based on an attractive U Hubbard model on a lattice with up to second neighbor hopping we 
derive an effective Hamiltonian for phase fluctuations. The superconducting gap is assumed to 
have s-wave symmetry. The effective Hamiltonian we finally arrive at is of the extended XY type. 
While it correctly reduces to a simple XY in the continuum limit, in the general case, it contains 
higher neighbor interaction in spin space. An important feature of our Hamiltonian is that it gives a 
much larger fluctuation region between the Berezinskii-Kosterlitz-Thouless transition temperature 
identified with T c for superconducting and the mean field transition temperature identified with the 
pseudogap temperature. 
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I. INTRODUCTION 

The origin of the pseudogap phenomena— in high-T c 
superconductors is one of the most challenging questions 
in the theory of superconductivity but there is no con- 
sensus on the correct theoretical approach to be taken 
to describe such phenomena. What is generally accepted 
is that the pseudogap is a manifestation of strong cor- 
relation effects which become progressively more impor- 
tant as the doping is reduced into the underdoped regime 
and the Hubbard Mott insulating state is approached. 
Roughly speaking, theories of the pseudogap state can 
be divided into two classes. 

One is based on the idea of a precursor state to super- 
conductivity. In this scenario the Cooper pairs are pre- 
formed below the pseudogap temperature (T*) assumed 
larger than the superconducting transition temperature 
(T c ) but there is no phase coherence in the temperature 
range between T c and T*m£ The phase fluctuations de- 
stroy phase coherence and consequently the superfluid 
stiffness is zero in this temperature regime. Therefore, 
this picture envisions that the two gaps ( superconduct- 
ing and pseudogap ) have the same origin. The precursor 
state idea is based on the results of angle resolved pho- 
toemission spectroscopy (ARPES) experiments, 4 which 
show a smooth evolution of the pseudogap into the super- 
conducting gap as temperature is reduced. This is also 
consistent with tunneling results 5 of Renner et al. who 
found a smooth evolution of tunneling characteristics as 
a function of temperature showing the gap filling but not 
closing as one goes through the superconducting transi- 
tion temperature. On the other hand data by Krasnov 
et al& on intrinsic mesa junctions show a separate super- 
conducting peak and a pseudogap hump. These data do 
not favor a pre-formed pair interpretation although su- 
perconducting and pseudogap both appear pinned to the 
Fermi energy. A different precursor scheme^ includes in 
the calculations finite momentum pairs which exist at any 
finite temperature below T* but do not contribute to the 
superfluid density made up of zero momentum Cooper 
pairs. Above T c there remain no zero momentum pairs 



and therefore no superfluid density. This approach con- 
stitutes a natural extension of BCS theory to the strong 
coupling regime and has many features that agree with 
experiments on the cuprates. For example pseudogap ef- 
fects become more prominent as the doping is reduced 
into the underdoped regime of the phase diagram. 

A second class of theories assumes that there is a sec- 
ond order parameter which competes with the super- 
conducting gap and that the microscopic origin of the 
pseudogap associated with this second order parameter 
is different from that of the superconducting gap. The 
competition between the two gaps manifests itself most 
importantly in the underdoped regime and so in this dop- 
ing regime the two gaps co-exist. As the doping is in- 
creased towards optimum and into the overdoped regime 
the second order parameter is weakened and could even 
vanish in the overdoped case. This feature is quite dif- 
ferent from the precursor scenario. A recent theoretical 
approach^ based on this line of thought has assumed that 
the pseudogap has a d-wave charge density wave order, 
the so-called d x 2_ y 2-density wave (DDW). An experimen- 
tally observed phase diagram— for the high-T c cuprates 
is consistent with this competing gap picture. However, 
a more recent theoretical workifi implies that the DDW 
model may not be consistent with tunneling experiments. 
Transport properties^! of the DDW state may also be 
used to test this picture. 

In this paper we deal only with the phase fluctuation 
picture. Cu02 planes play a crucial role in high-Tc super- 
conductivity and inter layer coupling is weak. Here it is 
assumed that the two-dimensional nature of the cuprates 
is an important feature of these compounds. For simplic- 
ity, however, we will consider only s-wave pairing on a 
lattice so that strictly speaking our work will not apply 
to high-Tc cuprates which are known to have ci- wave sym- 
metry. In this sense the present work is a first essential 
step in a realistic treatment of the cuprates. 

Noting that spontaneous breaking of continuous sym- 
metry is impossible in 2D^£ one may conceive 1 ^ that 
the superconducting transition, at least in underdoped 
cuprates, is not of the BCS kind but is rather related to 
the Berezinskii-Kosterlitz-Thouless (BKT) transition, 14 
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which is found in two-dimensional systems i^i^iL Of 
course, in reality, the coupling between Cu02 planes must 
become important right at T c since the transition must 
be three-dimensional with true off-diagonal long range 
order. 

The two-dimensional XY models on a square lattice is 
an example of the BKT transition. This model describes 
a system of spins Si constrained to rotate on the lattice, 
where i indicates the site. The Hamiltonian of the model 
can be written in the form —JJ2<ij> ^i'^ji where J > 
and < ij > means a sum over nearest neighbor pairs. It 
has been proposed^ that the BKT transition is associ- 
ated with unbound vortex-antivortex pairs. The mean- 
square distance (r 2 ) between a vortex and an antivortex 



<r 2 > = 



7T J/T - 1 

nJ/T-2 



(1) 



As temperature approaches the BKT transition temper- 
ature Tbkt — 7r J/2, (r 2 ) diverges and the pairs become 
unbound. 

In a theoretical approach to the precursor scenario, 
Cooper pairs are formed at the mean-field temperature 
(Tmf), which is identified as T* , The phase fluctuations 
destroy the long range order until the phases are locked 
in at the BKT transition temperature Tbkt, which is 
T c . Consequently, the effective Hamiltonian for phase 
fluctuations is mapped into a XY model Hamiltonian. In 
order to consider phase fluctuations of the superconduct- 
ing order parameter, one has to go beyond mean-field 
BCS theory. This can be implemented by using the well- 
established path-integral formalism for fermionic fields^ 
The effective phase Hamiltonian for an s-wave super- 
conductor can be described by Hxy = Jxy J dr(X76) 2 , 
where 9 represents phase fluctuations. To obtain a lat- 
tice version of the XY Hamiltonian one might discretize 
Hxy- However, it has been showni^ that discretization 
does not pick up all the important features of the phase 
Hamiltonian on a lattice, which will be derived and stud- 
ied in detail in this paper. 

The paper is organized as follows: In Section II, we 
introduce the path integral approach and derive an ef- 
fective action for phase fluctuations. We apply the 
Hubbard-Stratonovich transformation to integrate out 
the fermionic degrees of freedom. We also consider a 
case where the phase is treated in the continuum limit 
while the fermions still remain on a lattice. This serves 
as a bridge to the continuum limit. We show in section 
III that the effective Hamiltonian is not of the usual XY 



type. There are extra terms which do not belong to the 
usual XY Hamiltonian. It is shown how the new more 
complex Hamiltonian reduces to the XY model in the 
continuum limit. In this section we also explain how a 
periodic boundary condition applied to the lattice leads 
directly to a simple expression for the effective Hamilto- 
nian for phase fluctuations. In Section IV, we show that 
the effective Hamiltonian turns out to be of an extended 
XY type and investigate the effects of the extra terms. 
When we include next-nearest neighbor hopping in the 
electronic energy dispersion (Section V), we find that the 
extended nature of the Hamiltonian is preserved and the 
correction to the usual XY model becomes even more es- 
sential. We also explain a qualitative relation between 
the physics of the BCS-Bose Einstein (BE) crossovei—2. 
and the BKT transition temperature. Section VI con- 
tains discussions and conclusions. 



II. FORMALISM 

We begin with a 2D attractive Hubbard model on a 
square lattice in its simplest form conceivable, namely 



n 



<ij><7 



(2) 



where Ci„ is a fermion field with a spin a, /i is a 
chemical potential, U(> 0) is the pairing interaction, 
n-ia = Cf a Ci a , and t describes nearest neighbor hopping. 
To start we consider only the nearest neighbor hopping 
for algebraic simplicity; however, later we will include 
the next nearest neighbor hopping. The symbol < ij > 
means a sum carried out over nearest neighbor pairs as 
mentioned earlier. When we need a double summation 
over i and j, the usual notation ^ . will be used. We 
will set the lattice constant a = 1 and also use units such 
that h = ks = 1. The partition function (Z) of this 
model in the path-integral formalism is written as 



Z = J VC+VCcxp[-S] , 



(3) 



where the action S = J rfr{^ i(T C^d T Ci a + H} and the 
range for the integral over the imaginary time r in the 
action is from to 1/T. 

Defining spinors 'J',; for the fields Cjj and Cf, such that 
= (Ci Cii), the action S becomes 



S 




E 

<ij> 

r 



(4) 



where f a (a — 0,1,2,3) are Pauli matrices and f± = 



(fi ± if2)/2. The standard manipulations to derive the 
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effective action for phase fluctuations includes: i) the 
Hubbard-Stratonovich transformation^ to deal with the 
four-fermion-field term and ii) a gauge transformation 
to avoid violating the Mermin- Wagner theorem^ In or- 



der to use the Hubbard-Stratonovich transformation, one 
needs to introduce a complex auxiliary field <pi and the 
partition function Z becomes a function of ^Sf and <j>: 



where 



S = j dT \jjY. I^l 2 + E *i fo 8 * ^ *< - 1 E - E [<lH*t*+*< + M } 

K i i <ij> i J 



(5) 



(6) 



Note that if we require the equation for </>*; namely 6S[4>, 4>*]/S<p* = 0, the solution of this equation is <pi = UCiiC^. 
Consequently, 4>i can be identified as the order parameter. 

Parameterizing the auxiliary field 4>j — Aje l6j and making a gauge transformation tyj — exp[iT^8j /2}xi, where Xi 
is a spinor for neutral fermions, one can show that 



s = /H^EN a + E: 



r d T + T 3 ^(d T 9j) - nfi - TiAj 



E^i 



J,i+8 e 



-ir 3 6 itj /2 



(7) 



where O^j = 6i — 8j is a phase difference between sites i and j, S = ±x , iy for a square lattice. Note that if we equate 
4>j = UCjiCjj and <j>j = Aj-e*^' 2 , then we know that the phase part of Cj a is exp[i8j/2]. This means that the gauge 
transformation 21 - 22 splits the charged fermion field Cj a into a charge bose field exp[i0j/2] and a neutral fermion field 
Xj, which obeys a Grassmann algebra. 

After integrating out fermion fields, we obtain Z — J £>0*X>0exp[— S e ff], where the effective action is 



S eff = J dr^il&f-TrHG- 1 ] . 



(8) 



Here, Tr means the trace over the functional space, depending on the representation, and on the spin space. We will 
use a symbol tr for a trace over spin space only. In the real (lattice) space, the Green function can be represented as 



where 



G- 1 (i,j) = Gv 1 (i,j)-X(i,j) , 



5 



and the self energy has the form 



j) = --',u7.«, irt,. , + Utq ^ S 3,i+8 Sin ( ~lf ) + tf 3 E 5 3^+ 8 

5 



1 — COS 



'J_J_ 

2 



(9) 
(10) 

(11) 



r 



In the expression for Gq the bulk value of A is effects of phase fluctuations. For simplicity we also con- 
used. This means we consider only fluctuations of the sider only the static case (d T 9i = 0); therefore, we are not 
phase and assume no fluctuation in the magnitude of the concerned with the Landau damping effects which are re- 
gap. Magnitude fluctuation will also change the self en- lated to the scattering of the thermally excited quasipar- 
ergy £. However, in this paper we concentrate on the tides with the excitations of the phase field. However, if 
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one is interested in dynamics associated with the phase difference to the phase derivative as follows: 
fluctuations, then the Landau effects have to be taken 
into account^ 

Before going further we would like to see what hap- ^ 
pens if the phase is treated in the continuum limit while xj~7b& / J $j,i+6 s i n (^i,j/2)Xj — ^rp-it xt@i,i+8Xi+S 

the fermions are kept on a lattice. Let us consider jj s i,s 

jXi^{hj)Xj m the static case. Since we consider (12) 
the continuum limit for the phase, we change the phase Since 9i ; i + s = 9{ — 9 i+ s ~ —S ■ X79 — |(<5 • V) 2 9, 



~foitJ2Xi d is+sXi+s^~mY,Xi 6-V9+^{5-V) 2 9 X i+s ■ (13) 

i,8 i,8 

The fermions are on the lattice so we need to use Bloch state to expand Xi — e r< Xk« Now we have for example 

Y,*t(S ■ V9) Xl+ s = 2iJ2xi \{d x 9) sm(k x ) + (d y 9) sin{k y )] Xk . (14) 

i,S k 

Note that ^2$(S ■ V9) cos(k ■ S) — because of symmetry. We also have 

£tf" • V) 2 ^ +A - = [(d 2 J)cos(k x ) + (d 2 9)co S (k y )] Xk . (15) 

i,S k 



We therefore obtain 



Similarly, one can show that 



- f °X>k{ - \ lt M9)co S (k x ) + (dp)coa(k y )] 

- t [{d x e)$m{k x ) + (d y 6)sm(k y )] jxk . (16) 

H xt^t hi+d 1 - C0S (<W 2 )]xj 

~ xi [{d x 9f cos(fc x ) + (d y 9) 2 cosfo)] Xk . (17) 

k 

Consequently, when the phase is taken in the continuum limit while the fermions are still on the lattice, the self energy 
in the momentum space is 

£= f I —it [(d 2 x 9)cos(k x ) + (d 2 9)cos(k y )] +t[(d x 9)sm(k x ) + {d y 9) sm(k y )}\ 

+ \tf 3 [(d x 9) 2 cos(k x ) + (d y 9) 2 cos(k y )] , (18) 



which is the same as the self energy of Ref 



24 



III. EFFECTIVE ACTION tion used, the formalism associated with the derivation of 

the desired action has become quite simple because the 
In this section we concentrate on the effective action for self-energy that we have obtained in the previous section 
phase fluctuations. By virtue of the gauge transforma- depends only on the phase of the order parameter. 
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The effective action S e ff can be separated into the 
mean- field part (S^) and the phase fluctuation part (Sg) 
as follows: 

Seff = I drYthk? -TtHG- 1 ] 

i 

= 5 (0) -Trln[l-G E] 

= S (0) + Sg , (19) 

where 

S<°>= [drJ2^\ 2 -Trln[G^] , (20) 

i 

and 

oo 1 

5 e = TrV-[G S]" . (21) 
* — ' n 

n=l 

In the saddle-point approximation^ dS 1 - ^ /dA — re- 
duces to the BCS mean-field gap equation: 




and — (T/V)(dS^ ' /d/i), where V is a volume of the sys- 
tem, gives an equation for the filling factor: 

»=i-x;ft m h(§), (23) 

where E k = + A 2 with £ k = 

—2t [cos(k x ) + cos(ky)] — fx. In the calculation of 



and 



Since 



where r\ — > + is a convergence factor^ which origi- 
nates from G(0,t — r + ). Note that E^ ) makes no con- 
tribution to S^\ Since S^ include the phase difference 
between the nearest neighbor sites, the first non-trivial 



the effective phase-only action Sg, we assumed that 
the phase fluctuations are small so that it is sufficient 
to consider only the first and the second trace in the 
expansion for Sg; in other words, Sg ~ ffW + S&\ where 

S^ = Tr [G„E] (24) 

and 



S<- 2) = -Tr [G EG E] . (25) 



At finite temperature we obtain an extended XY model 
Hamiltonian which includes not only the nearest neigh- 
bor spin-spin interaction but also the next nearest neigh- 
bor and third neighbor interaction. This constitutes an 
extension of the usual XY Hamilton (Hxy), which con- 
tains only nearest neighbor interactions. Inclusion of 
the next nearest neighbor hopping in the electron energy 
dispersion makes the extended feature of the effective 
Hamiltonian even more robust. In this case the effec- 
tive Hamiltonian manifests this extended property even 
at zero temperature. This will be discussed in section V. 

The easiest way to calculate Sg is in the momentum 
representation for Go and E. For simplicity, we will use 
a four vector notation: K = (k, w n ), J2k = Sk> 
where oj n is a fermionic Matsubara frequency, and 
J dXi = J dr^i with the sum over the lattice sites. 
It can be shown that (K\G\K') = S(K - K')G(K') and 
(K'\E\K) = t &°\K' -K,\s.)+f 3 & 3 \K' -K,k), where 



(26) 



(27) 



(28) 



term is 8fj. We obtain immediately the action of the 
usual XY model; however, a Of ^ term will also appear 
in S^ 2 \ Moreover, the phase difference between the next 
nearest neighbor sites is also obtained from the second 



&°\K' -K,k)=it J dX t e 



J^e*-' sin 



U,i+S 



( K < - K ,k) = t f dXie-^'-V** ' 
teT,K K>( K \ G o\ K '}( K 'm K )' we obtain 



ik-6 



1 — COS 



'i,i+S 



K 





\- cos ( e jA~ 

V 2 )_ 


<ij> 
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trace. As we will see later Section V, if we include the The second trace can also be calculated in a similar 
next nearest neighbor hopping, then even contains manner even though it is more complicated: 
the phase difference between next nearest neighbor sites. 



S {2) = -tr G{K)t{K - K', k')G{K')E(K' - K, k) , (29) 

K,K' 

where S = foS^ + f 3 & 3 \ After some manipulation, one can arrive at 

S(2) = ^ E tr i G ( K ) G ( K '^ S (0) (K - K\ k')£ (0) {K 1 - K, k) 

K,K' 

+ \ tr[G(K)r 3 G(K')r 3 }t^(K~K',k')t^(K' -K,k). (30) 

K,K' 

It can be shown that the mixed terms do not contribute to in our consideration which deals only with a local 
effective action. In order to calculate we need further manipulations. Define A$(Xi) — sm(di t i+s /2) and express 
it in terms of its momentum counterpart as follows: As{Xi) = J2q As(Q) expfiQX;], where Q is a four vector (q, n m ) 
with a bosonic Matsubara frequency f2 m . Then, for example, one can show for the first term of that 

J2 tr [G{K)G{K')] £ (Q \K - K' , k')t (0 \K' - K, k) 

K,K' 

= (tt) 2 E E ^(«^f Q)^ ,q "' E tr [G(Jf )G(K - Q)] e*^^') . (31) 

5,5' Q K 

I 



See Appendix A for a detailed derivation. A similar ex- 
pression can be obtained for the second term. Since, as 
we mentioned earlier, we are concerned only with the lo- 
cal effective action for the phase only and since we wish 
to see if the equivalence between the phase-only Hamil- 
tonian and the usual Hxy, which is found in the contin- 
uum limit, still holds on a lattice, we expand about 
Q = and keep the leading order. Here we would like 
to explain that in contrast with the continuum limit this 
procedure requires a periodic lattice. In the continuum 
limit As(Q)As'(— Q) ~ q 2 8 q 9- q ; therefore, we can put 
Q = for tr [G(K)G(K - Q)}. However, this is not the 
case on the lattice because Ag(Xi) does not deal with the 



phase derivative but rather it deals with the phase differ- 
ence. One might expand A$(Q) and tr [G(K)G(K — Q)] 
in terms of Q and, then, regroup terms of the same order 
of Q. If we do so however, we cannot transform back to 
the lattice space with a function of the phase difference 
Instead of doing so, we see under what condition we may 
simply put Q = for tr [G(K)G(K ~ Q)] and neglect 
corrections. What we found is a periodic boundary con- 
dition for the phase 0i. If the size of the 2D lattice is 
L x L, then 9(jCi + L) = 8(ri). See Appendix B for de- 
tails. 

Now, it can be shown that 



1 5,5' K J 



, 'i,i+5 \ ■ / 

sin sin 



V E E tr [G(K)r 3 G(K)f 3 ] e ik 'CS-0 



5,5' K 



X / dX 4 



1 — COS 



'i,i+5 



1 — COS 



'i,i+5' 



(32) 



In order to compare the effective Hamiltonian with the Hxy, we need to expand in terms of the phase 
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differences whiie keeping terms up to 0? .. Since the 
first non-trivial term for the second term is 6fj, we 
will ignore it in the effective Hamiltonian. Using sym- 
metry properties for G(K) such as remaining the same 
with respect to the exchange of k x and k y , we obtain 
the effective local Hamiltonian: He = (6\M\6), where 

(9\ = (0i : i+x,6i,i~x,6i,i+y,0i,i-y) and 



M 



( a /3 7 7 

(3 a 7 7 

7 7 a f3 

V 7 7 /3 a 



with components as follow: 



a = 1* £ tr [G(K)f 3 e^3] cos(k x ) + \t 2 £ tr [G(K)G(K)} 



A 



P = it 2 ^tr[G(/OG(/0]cos(2fc x ) , 



K 



A" 



(33) 
(34) 



and 



l=^t 2 J2 tr [G(if)G(X)] cos(fc a ) cos(fc y ) 



(35) 



A 



As one can see, the 4x4 matrix M is not diagonal. This 
means that the effective Hamiltonian He is not equivalent 
to the usual Hxy, which would be pure-diagonal. Not 
only do we have terms like 0? i+ ~ and 0?j_ s but we also 
have terms of the form See Appendix A. We 

will investigate effects of these off-diagonal terms in the 
next section. 

It is worthwhile making sure our new result for He 
reduces to the well known XY-type Hamiltonian^^ in the 
continuum limit. For this purpose, we need to recover 
the explicit dependence on the lattice constant a in the 
expression for He which we have set equal to 1 . This will 
allow us to track orders of a and take the limit a — ► 0. In 
this case the phase difference becomes a derivative of the 
phase (V0) while ta 2 — ► where m is an effective mass 
of the electron. It is straightforward to show that, in the 



continuum limit with £k 
where 



2m 



H, He 



H, 



(i) 



H, 



(2) 



(36) 



and 



H. 



(2) 



16to 2 



J2^[G(K)G{K)]k 2 / dr(V6l) 2 . (37) 



A 



See Appendix C for details. Consequently, the effective 
Hamiltonian He we have derived does reduce to the well- 
known XY-type expression, H X y = Jxy J dr (V6>) 2 . 

IV. OFF-DIAGONAL TERMS 

Let us now consider effects of the off-diagonal terms in 
He- At T = 0, He again becomes equivalent to the XY- 
type Hamiltonian because the off-diagonal terms of M, 
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FIG. 1: Geometrical illustrations of induced interactions for 
spin pairs. 5 is a vector to the nearest neighbor site, 82 is to 
the next nearest neighbor site and so forth. 



namely (3 and 7 vanish. However, at finite temperature 
they are finite so that He is no longer of the usual XY 
type; that is 



He — a ^ 9 2 j + 2/3^ (Qi,i+x0i,i-x + 6i,i+yQi,i-y) 
<ij> i 

+ 27 [6j t j +x + 6j t j-x) {Qi,i+y + @i,i-y) ■ (38) 



Since the phase fluctuation between two sites is small 
by assumption, one can make the approximation that 
di^i+x — sin^j-i-^). Introducing a 2D classical spin 
Si = (cos(0j), sin(#i)) at a site i, it can be shown, 
within the approximation we have made, that, for ex- 
ample, t ,i+x0i,i-x — (Si x S l+£ ) ■ (Si x Si_*). It can 
be seen that a spin at the site i + x couples to a spin 
at i — x, which is the next nearest neighbor spin-spin in- 
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teraction. The only assumption we make is that terms 
higher than 0? . for a given are negligible. This 

does not mean, however, that only the nearest neighbor 
phase differences are important. Following the proce- 
dure we described in Appendix D, one can show that the 
second term of Tig, which is proportional to /3, becomes 
P H<ij> Qlj+I 3 J2<i J>3 S i' S j U P t0 a constant, where the 
symbol < ij >s indicates a sum over the next-next near- 
est pairs. Similarly, the third term of Jig proportional 
to 7 turns out to be 2jY,<ij > 9 i,j + 2 7E« y » s i ' s j> 
where -C ij 3> means a sum over the next nearest pairs. 
Since 0? ■ can also be represented in term of Si • Sj, the 
effective Hamiltonian Tig can be written as 

Tig = —Ji ^2 ' &i + ^ 2 E ^ ' Sj + ^3 E ^ ' ' 

<ij> <S»i» <U>3 

(39) 

where J x = 2(a + /3 + 27), J 2 = 27, and J 3 = /3. It is 
clear that this Hamiltonian is not of the usual XY type 
but instead of an extended XY type. A geometrical ex- 
planation for the appearance of the next nearest neighbor 
and the next-next nearest neighbor term in Eq. I|39l) is 
illustrated in Fig. 1(a). Indeed, these terms come from 
the second trace proportional to t 2 , which has a factor 
E(°)S(°). Each self-energy picks up <5 = ±x, ±y, 
and the second trace gives terms with resulting vectors 
8 + 8' = 62 or S3 , where 82 — ±x ± y and £3 ± 2i, ±2y. 
This geometrical picture also works when we include the 
next nearest neighbor hopping (£') in the electron disper- 
sion curves. 

The physics of Tig depends on the relative magnitudes 
of coefficients Ji, J2, and J3 in Eq. l|39fl as well as their 
relative signs. For example, if f3 is negligible, and a and 
7 are both positive so that J\ > 2J2, Tig describe a non- 
frustrated XY model, and its critical behavior can be un- 
derstood in terms of the usual XY Hamiltonian with an 
effective coupling constant J e ff = (J\ — 2 J2)~ However, 
this does not mean that in this case Jig is equivalent to 
~Jeff ^2<ij> Si 1 Sj because the local behavior of these 
two Hamiltonians is different. In general, however, as 
long as J\ is dominant, the large length scale behavior 
is of the usual XY type. To calculate Ji, J2, and J 3 , 
we need to know how A and fj, change with increasing 



temperature. We choose the pairing interaction U = lAt 
and the filling factor n — 0.9 and self-consistently solve 
Eqs. (E2Jl and (O to determine A(T) and fi(T) as func- 
tions of T. In Fig. 2 we plot J e ff vs T scaled by Tmf 
for t' = 0. The value of the BKT transition tempera- 
ture (Tbkt) is indicated by an arrow. As one can see, 
in this case Tbkt/Tmf — 0.73. If n gets closer to 1 as 
well as U becomes smaller; namely, n ps 1 and U < t, 
Tbkt/Tmf ~ 1 as we verified numerically. 



2.0 




FIG. 2: T^Jeff (solid curve) as a function of temperature (T) 
scaled by Tmf for the pairing interaction U = lAt and the 
filling factor n = 0.9 with next nearest neighbor hopping t! = 
0. The solution of T = %J e ff(T) gives the BKT temperature 
Tbkt which is indicated as an arrow. 



V. INCLUSION OF THE NEXT-NEAREST 
HOPPING 

When we include the next nearest neighbor hopping 
(t'), the Hamiltonian becomes 



n = t c+c jo -a c&c ia M E -t/£ n, 



(40) 



<ij>(7 



<ij No- 



where <C ij 3> denotes a summation over the next nearest neighbor sites, and now £k = — 2t [cos(k x ) + cos(k y ) 
At' cos(k x ) cos(ky) — fi. The same procedure as before yields the actions 



T <9 r + T 3 -(d T 9j) - f 3 ^i - fiAj 



s i,jXj 



-tf 3 S 3 



-if 3 e i , j /2 



E 5 3 



i+s 2 e 



-if 3 9 itj /2 



Xj 



(41) 
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and the Green function 



i+s 2 

. s s 2 



with the self energies 



(42) 



'—^itJ2hi+s* sin' 



(43) 



and 



& 3 \i,j)=t 53 <w 



1 — cos I 

2 



1 — COS I 

2 



(44) 



We obtain the effective Hamiltonian including the next nearest neighbor hopping contribution: 

T~Le = —Ji 53 • Sj + J2 53 S« ' ^ ^ 3 5Z ' ^ 4 ^' ' ^ 5 ' ' 

<ij> -Cij> <ij>3 <ij>i <ij>5 



(45) 



where V^w,-^ means V, A . Figs. 1(b) and (c) show 

■^<y>4(5) ^— 'M4(5) ° V ' V ' 

geometrical descriptions for theses terms. In the same 
way that the t 2 term gives interactions of the form S; ■ 
Si+s 2 and 8^ • Si+5 3 , the tt' term induces the form Si ■ Si+5 4 
because S+S2 = 64, and the t' 2 term gives rise to • S^+^g 
and Si ■ Si + g 5 because S2 + S' 2 = S3 or #5, where 64 — 
±x ± 2y, ±2x ± y, and 5$ = ±2x ± 2y. 



It is obvious that the extended features of Tie are ro- 
bust in band structure parameters. Moreover, contrary 
to what we found before, now J2 is finite even at T = 
because the i'-term contributes to the first trace as well 
as to the second trace as shown in Appendix E, where 
detailed expressions for Ji (i — 1, 2, 3, 4, 5) are presented. 
In order to calculate the J7i's we choose t' = — 0.2£ for sev- 
eral values of parameters ([/, n). In Fig. 3, from the top to 
bottom panel (U, n) are (1.4, 0.9), (2.0, 0.9) and 2.0, 0.4), 
respectively. Since our numerical calculation indicates 
that J\ and J2 are dominant and J\ > 2J2 in most of 
the temperature range, we neglect J3, J4, and J5. In this 
case the critical behavior is still of the usual XY type so 
that we introduce an effective coupling J e ff = J\ — 2J 2 . 
As also shown in Fig. 3, (it /2)J e f f /Tmf becomes smaller 
as U is increased or n is decreased; therefore, the cross- 
ing point between (it /2)J e ff /Tmf and T/Tmf changes 
correspondingly. For the attractive Hubbard model, the 
(U,n) phase diagram shows BCS, cross-over, and BE 
regimes; for example, a large n(« 1) and a small ?7(<C t) 
corresponds to the BCS regime while a small n and a 
large U corresponds to the BE regime^ Consequently, 
in a qualitative sense Tbkt/Tmf becomes smaller as the 
parameter moves from the BCS to the BE regime. 



VI. DISCUSSIONS AND CONCLUSIONS 

An effective Hamiltonian for phase fluctuations has 
been derived on a square lattice using the attractive Hub- 
bard model to describe electron dynamics. We find that 
the effective Hamiltonian is not of the usual XY type, 
which maps onto a spin Hamiltonian with nearest neigh- 
bors only, but is extended XY containing spin interaction 
between second up to fifth neighbors. This is in contrast 
to the common assumption that discretization of the con- 
tinuum limit Hamiltonian leads to its lattice version. In 
the continuum limit the BCS phase stiffness and hence 
the effective Hamiltonian vanish at Tmf because they 
are proportional to the number of Cooper pairs. In other 
words Tmf can be determined not only from the mean 
field equation for the gap but also from the phase stiff- 
ness in the continuum limit. However, this is not the case 
on the lattice. Since we showed that the effective Hamil- 
tonian we obtained reduces properly to the well-known 
continuum limit result, this is clearly a special property 
of the continuum limit and does not hold on a lattice. It 
is interesting to note that the periodic boundary condi- 
tions in our formalism are essential to get simple results 
and they are those that are usually assumed for summa- 
tion over the Brillouin zone. Without them the effective 
Hamiltonian would not necessarily end at a finite number 
of neighbors. 

The extended feature of the effective phase Hamilto- 
nian is reinforced when the next nearest neighbor hop- 
ping is considered in the electronic dispersion. The crit- 
ical behavior is however still BKT-like but with an ef- 
fective coupling which depends significantly on temper- 
ature. This arises because near the BKT transition the 
extended Hamiltonian reduces to the XY case but this 
simplification only holds in this restricted region of tem- 
perature. In particular the generalized phase Hamilto- 
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FIG. 3: §<Jeff (solid curves) as functions of temperature 
for various values of the pairing interaction U and the fill- 
ing factor n. From the top to bottom (U,n) are (1.4,0.9), 
(2.0,0.9) and 2.0,0.4), respectively, with t' = 0.2t. The BKT 
temperature Tbkt is indicated as an arrow in each panel It 
is given by the intersection of the solid curve for \J & ff and 
the dashed line, (ir /2)Sf e f / /Tmf becomes smaller as U is 
increased or n is decreased; therefore, the crossing point be- 
tween (n/2)J' e ff /Tmf and T/Tmf changes correspondingly. 
As we move away from weak BCS regimes toward strong cou- 
pling the size of the phase fluctuation regime between the 
mean field temperature (T*) and the BKT temperature (T c ) 
increases. 



nian does not vanish at the mean field temperature Tmf 
which is identified with the pseudogap temperature T*. 
The crossover from BCS to Bose-Einstein corresponding 
to increasing U and/or decreasing n reduces the ratio of 
the BKT transition temperature to the mean field pseu- 
dogap temperature. 

In the continuum limit, 1 — Tbkt /Tmf — 4Tmf/ Ef, 
where Ep is the Fermi energy, for the BCS limit (large 
n and small U)~ Since Tmf/Ef is negligibly small for 
typical values of T MF » 20-30 meVand£ F « 800-1000 
meV, the temperature region between Tbkt{= T c ) and 
Tmf{= T*) is negligibly small as well. However, on the 
lattice we find instead 1—Tbkt/Tmf — 0.375 for n = 0.9 
and U = 1.4, which are typical values of these parameters 
in the BCS regime. Beyond the BCS regime (small n and 
large U), Tbkt /Tmf is mostly 0.5 in the continuum limit 
while the ratio reduces to 0.22 for n = 0.4 and U — 2.8 on 
the lattice. Therefore, the range for phase fluctuations 
between T c and T* can be much larger on the lattice as 
compared with the continuum limit. 

In our consideration J\ (or J\ for t' ^ 0) is not the 
BCS phase stiffness. Thus it does not have to vanish at 
Tmf- However, J e // (or Jeff) is the phase stiffness near 
the critical region and it determines Tbkt- Nevertheless 



we want to mention that we cannot interpret J e f / as the 
BCS stiffness in the whole temperature range. It is only 
near Tbkt that ~JeffJ2<ij> &i ' &j can represent the 
full Hamiltonian Tig that we obtained on the lattice. 

We have considered only s-wave pairing in this paper. 
For cuprates d-wave pairing is appropriate and this case 
received considerable attention in the literature but only 
in the continuum limit ^S^S^S^S^L In the continuum 
limit the order parameter </>(R, r), where R(r) is the 
center of mass (relative) coordinate, is represented by 
</>(R, r) = A(R, r)e* e ( R - ) because the bond phase is usu- 
ally replaced by the site phase. However, on the lattice 
it is the bonding order between the nearest sites that 
we must use as the building block to construct the d- 
wavc order parameter. This means that we can consider 
phase fluctuations between the nearest neighbor building 
blocks in addition to the phase difference by ir between 
the two; therefore, we may expect that a subdominant 
component to the order parameter may be induced by 
phase fluctuations. 

We neglected the so-called Landau terms associated 
with damping effects in this work. These terms do not 
enter in the static case which we have considered here. 
For dynamics, the Landau damping effects may play a 
role. Nevertheless these effects have been shownSi not 
to be important for T < O.GTmf in the case of s-wave 
superconductor in the continuum limit. However, it has 
also been pointed out that because of the nodal struc- 
ture of the order parameter of a d-wave superconduc- 
tor the Landau terms can have important effects even 
at low temperature i24 Therefore, to describe the dynam- 
ics of the phase fluctuations in a d-wave superconductor 
one has to consider Landau damping effects. Such ef- 
fects will also be important on a lattice. An extension 
of the present work to include Landau terms is necessary 
if dynamics is considered. Further complications associ- 
ated with the existence of nodal quasiparticles have been 
noted in Refi^i 
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APPENDIX A: MANIPULATIONS FOR THE 
SECOND TRACE 

In this Appendix we manipulate the second trace. Let 
us consider 
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53 tr [G{K)G{K')\ E (0) (K - K', k')£ (0) - K, k) 



K.K' 



where 



with 



= sin 



Y,MQV QXi 



(Al) 



(A2) 



(A3) 



Now Eq. IjA 111 become 



53 tr[G(/f)G(K")] / dXidXje-^-K'XX'-X') 

K,K' J 

x(it) 2 Y,e*'- s+ik - 5 ' ^ A s (QyQ x >A s ,(Q'yQ' x i 

8,8' Q,Q> 

(it) 2 53 53 MQ)A 5 , {-Q)e~^ s 53 tr [G(K)G(K - Q)] e < k ^+ 5 ') 



(A4) 



8,8' 



K 



For the local effective action, we put Q = in tr [G(K)G(K — Q)], which means that we keep only the leading 
term in the expansion of Eq. IjAlj) . 1 '^ 2 ?' 24 Note that sin^i+a) ~ 9^i+g ~ V# ~ qO(q) so that sm(6i^ + g) sin^j^+j/ ) ~ 
q 2 6{q)Q(—q). Consequently, tr [G{K)G{K — Q)] with a finite Q gives higher order terms in the expansion. Since 



Y,MQ)A&'{-Q)e-^- 5 = J dX t A s (t-6,T)A s ,(i,T) 



we have for Eq. (|A1(I 



t 2 ^^tY[G{K)G{K)]^< s - 5 ") f dX t 

s r/ is- J 



8,8' K 



, 'i,i+S \ ■ / Vi.i+S' 

sin — - — sm 



(A5) 



(A6) 



The manipulation for the second term exactly parallels the above derivation. However, as we can see, the second term 
is negligible because of the assumption of slow variation of the phase; therefore, we neglect it. 
Now we have an effective Hamiltonian: 



We// = tY,^[G(K)f 3 e l ^ f3 ]cos(k x )T 53 

K <ij> 



1 — cos 



%3 



8,8' K 



sm — ■ sm 1 



(A7) 



For a small variation of one can derive an effective local Hamiltonian Tig. Its derivation is straight- 

forward so that we only briefly mention some details about it. For the second term of Tt e ff, if (5 = 8', 
then one has itr [G(-ftT) 2 ] J2i s (^« ~~ @i+s) 2 - When 6 / 6', for example 8 — x and 6' — —x, one has 
jtv [G(K) 2 ~\ cos(2fc x ) Y^,i {Si — @i+x) (#j — 6i-x)- For other cases, one can apply similar manipulation. Note that 5 
and 5' are exchangeable. 



APPENDIX B: BOUNDARY CONDITION 



Let us explain how we end up with the boundary condition. Consider 

.0) 

,K-Q 



53 53 mq)a s , (-Q)e-^ 53 e *W) Aj$ 

8,8' Q 



K 
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A (00) I 2^ A (00) 

A w + 2 Q dQ 2 



Q=0 



E E Mq)a*> i-Q)^- 5 E ^ (W,) aS 

5,5' Q K 



+ EE^w^-^^E 6 

5,5' Q 



K 



2^ 



5Q ; 



.a(°°) 



+ 



(Bl) 



Q=0 



where A^°^_g = tr [G(fc)G(i ; ir — Q)}. In the bracket [ ], we do not have odd power terms in Q because of symmetry. 

The first term contributes and gives 0ij + s9 i+ s' to S^ 2 \ Now, we need to see under what condition the remaining 
terms are negligible. Let us transform back to the lattice space using 



MQ) = J dXiAsiXJe-^ 



(B2) 



Then we have, for the second term, 



E E MQ)** (-Q)e-^ s E ^ {s+s,) Iq 2 



5,8' Q 



K 



dQ 



-A(oo) 



Q=0 



= f dXidxl E E MA)A S > {vi + sy k 

' 5,S> Q,K 



(6+S') 1 d 2 jQlXt-X',) 

2 dXidX' 



dQ 



_A(oo) 



(B3) 



Q=0 



If the lattice size is large enough that we can change ^\ ~* I dri, then assuming a periodic boundary condition 
9{vi + L) — 6(ri), we have 



J dXidXlAsir'JAs'iri+S) 



d 2 



dXidx' 



iQ(Xi-X<) 



—Asir^Ai^Ti + S) 



boundary 



In general if 



d 
dX~. 



As{ri)As>{ri + 6) 



boundary 



= o 



(B4) 



(B5) 



we may simply put Q = in A^°^-_q for S^ 2 K Since 9{vi) is periodic, so are its derivatives. Consequently, the above 
condition is satisfied by the periodic boundary condition. 



APPENDIX C: REDUCTION TO THE CONTINUUM LIMIT 



In this appendix we show how Tig reduces to the usual form in the continuum limit. In order to do so, we need 
to recover the lattice Constance a and take the limit a — > while keeping ta 2 — ► where m is the electron mass, 
because £ — > k 2 /2m — \i. Let us consider first 



H 



(i) 



Now consider H^: 

H (2) 



1 — COS 



'■'■■J 
2 



*E tr [G(K)f 3 e lr,ul " f3 ] cos(k x a)T E 

K <ij> 

t E tr [G(K)f 3 e^] (l i(k x a) 2 • • •) ± £ & ~ <W 
Etr[G(i^)f 3 e lT ' w " f3 ] f 

IS J 



8m 



(CI) 



>EE tr W 2 ] e w,) E 



sin 



8,6> K 



H,i+6 



i+b' 



- k E tr [OT) 2 ] £ (l - y [k • (5 - S')f\ £ (ft 

K S,S' ^ ) i 

- iE tr ^) 2 ] [dv[kl{d x 6f + kl{d y 6f] 

if J 



Oi+s) (Oi — Oi+s>) 



16m 2 



APPENDIX D: OFF-DIAGONAL TERMS 



The effective Hamiltonian is 



He — a £ Q 2 - + 2(1 £ {Ois+ x Qi,i-x + 9i,i+y9i,i-y) 

<ij> i 
+ {6j,i+x + 6j,i-x) {Qi,i+y + Qi,i-y) ■ 



First let us consider the second term. 

2/3 ^ ^ {Qi,i+x@i,i — x ~l~ @i,i-\-y y ) 

- 2 /3 £ [ sin (^,»+*) sin^^j-s ) + sin(0 iii+s ) sin(6> i;i _ s )] 

~ 2/3 £ [(Si x S<+*) • (S< x S<_a) + (S< x S i+S ) • {Si x Sj_ fl )] 

Applying a vector identity to this expression, we then obtain 

(Si x Si+j) • (Si x Si—x) = Si+x • Si—x ~ (Si • Si+x)(Si • Si—f) , 
where we used the fact that Sj ■ Sj = 1, we know that the second term becomes 

(3 £ Sj ■ Sj - 2/3 [cos(6 iti+x ) cos(6 iti - x ) + cos((9^ + g) cos(ft,»_g)] 

<ij>3 i 

Since cos(^i^ + i) ~ 1 — (0i^ +x ) 2 /2, we obtain, for the second term, 



2/3^ 



H^i+xvi^i — x 



<ij> <ij>3 



It is also straightforward to show that 



APPENDIX E: COUPLING CONSTANTS 



Ji = j{tJ2 tr[G(K)T 3 e">^ f3 ] cos(k x ) + t 2 J2 tr[G(K)G(K)] 

't' 



K 



K 



1 + cos(2fc x ) 



+ 2 cos(fc x ) cos(k y ) + 2 ( — ) cos(2fc x ) cos(2A: y ) 



14 



where rj — > + is a convergent factor. 

^ 2 = II -t'^te[G{K)Ti]cos{k x )cos{k y ) 

{ K 



K 



cos(k x ) cos(fcy) — 2 ( — J (cos(fca;) + cos(2/ca;) cos(fc y )) 



+ [ j) (1 + cos(2fc x ) cos(2fc y ) + 2 cos(2fc x )) 



(E2) 



k L ^ ' 

J 4 = jtt'^2tT[G{K)G(K)]cos(2k x )cos(k y ) 



cos(2fc x ) 



and 



In actual calculations, we use 



and 



J 5 = ii ,2 ^tr[G(K)G(K)]cos(2fc ;E )cos(2^) . 



A' 



tr[G(k, w n )f 3 e"""" f3 ] = !-■&■ tanh ' U 



£ k V 2T 



T^tr[G(k,c n ) 2 ] = -i 



2T 

r 



1 - tanh 2 



2 ( £k 

2T 



(E3) 
(E4) 

(E5) 
(E6) 
(E7) 
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